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of the circles M and N, go into A' and B', intersections of M a diameter of 
K' , and N' a circle diametric to K' . Hence A' and B' are inverse with re- 
spect to K' . 

This answers question (2). And question (1) may then be answered 
in the affirmative in view of the two following facts: 

First, any arbitrary finite group of fractional linear substitutions may 
be obtained by transforming one of Klein's canonical groups by a fractional 
linear substitution; and 

Secondly, if the poles of a fractional linear substitution are skew-in- 
verse with respect to some circle C, then this substitution corresponds, by 
stereographic projection, to a rotation of the sphere which has C as a great 
circle. 



HISTORICAL NOTE ON THE NEWTON-RAPHSON METHOD 
OF APPROXIMATION. 



By FLORIAN CAJORI, Colorado College. 



Newton explained his method of approximation to the real roots of 
numerical equations in a tract, De arwlysi per aequationes numero terminor- 
aum infinitas* which is celebrated chiefly as containing the first announce- 
ment of the principle of fluxions and the binomial theorem. Newton placed 
it in the hands of his teacher, Isaac Barrow, in 1669. Barrow sent it to John 
Collins, a member of the Royal Society, who, for his zeal for collecting and 
diffusing scientific information, received the sobriquet of "English Mer- 
senne." The tract became known to some correspondents of Collins and to 
friends of Newton, but it was not printed until 1704 and 1711. Substantial- 
ly the same explanation of his method of approximation to the roots of nu- 
merical equations was given by Newton in a second tract, the Methodus flux- 
ionum et serierum infinitarum, which was planned for publication in 1671, 
but was not printed until 1736. 

The earliest printed account of Newton's method of approximation 
appeared in Wallis' Algebra, London, 1685, chapter 94. Wallis explains 
Newton's process of solving y % — 2y-5=0, and derives, in the equation 
y a +axy+aay— x 3 — 2a s =0, Newton's value of y expressed in terms of a and 
a; by a rapidly converging series. The two equations are treated by Newton 
in much the same manner. The solution of y* — 2y— 5=0 is exhibited in 
Wallis' Algebra, as also in the two tracts of Newton referred to above, in 
paradigms which are identical except in some of the last digits in the deci- 
mal fractions. The following is copied from page 268 of Vol. I of Newton's 
Opera (Ed. Horsley) : 



* Isaac Newton's Opera (Edition by Horsley), 1779-1785, Vol. I, pp. i 
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y 3 -2y-5 


-0 


+2.10000000 
-0.00544853 


+2.09455147=y 


2-fp=y 


+v 3 

-2y 
-5 


+8+l2p+6p*+p* 

—4-2p 

-5 


Sumraa 


-l+10p+6p 2 +p :i 


0,l+q=p 


+F H 
+6p 2 
+10p 
-1 


+0,001+0,03g+0,3<f +g 3 
+0,06 +1,2 +6,0 
+1 +10, 
-1, 


Summa 


+0, 061+11, 23q + 6,3q 2 +q 3 


-0,0054+r=g 


+6,3g s 

+ll,23g 

+0,061 


+0, 000183708-0, 06804r+6,3r'" : 

-0,060642 +11,23 

+0,061 


Summa 


+0,000541708+ll,16196r+6,3r 2 


-0,00004854+s=r 







We see that Newton takes as the initial approximation to the real root 
y =---2. He expresses himself thus: "Sit 2 numerus qui minus quam decima 
sui parte differt a radice quaesita. Turn pono 2+p=y..." The equation in 
p becomes p a +6p 2 +10p— 1=0. Neglecting the higher powers of p, he gets 
10p-l=0; taking p=.l+q, he gets g 8 +6.3g 2 +11.23g+. 061=0. From 
11. 23g+. 061=0 he obtains q=— .0054+if, and by the same process, 
r = - . 00004853. Finally, y=2 +. 1 - . 0054 - . 00004853=2. 09455147. 

Newton seems quite aware that his method of approximation may fail. 
If there is doubt, he says, whether p=.l is sufficiently close to the truth, 
find p from 6p 8 +10p— 1=0, but he does not inquire whether even this latter 
method will always answer. He gives an explanation of the solution 
of y % —2y—5=0, but does not develop his method further. No other exam- 
ples are given. 

Wallis does not praise Newton's method over the older; he merely 
states that it "is very different from that of Vieta, Oughtred, and Harriot, 
which is commonly received." The words "very different" are in marked 
contrast to the statements of the modern historians, H. Hankel* and 
M. Cantor, f who make Vieta' s method of approximation appear as almost 
identical with the procedure given by Newton. The two are not the same. 
The essential difference lies in the divisors used in finding the successive 
approximations. If r is the approximation already reached, then Newton 
uses a divisor which in our modern notation takes the form /' (r) ; Vieta's 
divisor may be expressed in the form |/(r+s,.)— f(r) | — Si", where f(x) is 
the left side of the equation f(x) =k, n the degree of the equation; and s, is 



* H. Hankel, Geschichte der Mathematik, Leipzig, 1874, pp. 369, 370. 

t M. Cantor, Vorlesungen iiber Geschichte der Mathematik, Bd. II, Leipzig, 1900, pp. 640, 641. 
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a unit of the denomination of the digit next to be found. Thus, in Vieta's 
example, x*—5x*+ 500* =7905504, if r=20, Vieta's divisor is 878295, while 
Newton's divisor is 794500. The inaccuracy of Hankel and Cantor consists 
in attributing to Vieta the same divisor as that used by Newton. This error 
probably arose from the fact that in some quadratic equations given by Vieta 
the two divisors happen to be the same in value. The advantage of 
the Newtonian divisor over that of Vieta lies in the smaller amount of com- 
putation usually called for. 

The study of Newton's solution of y s — 2y — 5=0 reveals the fact that 
his method of approximation is not the same as what is called "Newton's 
method" in modern text-books. If r is the initial approximation to a root 

fir) fir ) 

of the equation f(x) =0, then r 1 --r— „ , \ and r.>=r, — i,\ x \ , are desig- 

f (r) f (rj 

nated in modern manuals as approximation by "Newton's method." But 
from the standpoint of the practical computor, this procedure is quite differ- 
ent from that of Newton. Newton derives, as we have seen, each succes- 
sive step p, q, r of approach to the root, from a new equation; while in the 
modern process just described each approach is made by substitution in the 
original equation. 

This modification of Newton's process was first effected by Joseph 
Raphson (1648 (?)— 1715(?)), a Fellow of the Royal Society, who published 
in 1690 in London, a booklet, bearing the title, Analysis aequationum 
universalis. A. second edition, with an appendix, appeared in 1697. * 

DeMorgan remarks that Raphson does not mention Newton ; Raphson 
evidently considered the difference sufficient for his method to be classed in- 
dependently. He does not mention Vieta. In our opinion, Raphson's paper 
is of much greater historical interest than is ordinarily attached to it. Its 
importance lies in the modification which it makes of Newton's method. 
Lagrange recognized the modified process as ' 'plus simple que celle de New- 
ton;" and it has now supplanted Newton's method. 

In view of the facts it is doubtful whether the method of approxima- 
tion described by Raphson should be named after Newton alone. In the 
first place, the prcesses used by Newton, though not identical to that of 
Vieta, resembles it. Newton merely simplified the divisor used. Hence the 
honor of invention falls largely on Vieta. In the second place, Newton did 
not develop his method further than simply to solve the cubic y % — 2y— 5=0. 
That Raphson worked independently of Newton we doubt. But Raphson's 
version of the process has been accepted as an improvement. It would seem, 
therefore, that the "Newton-Raphson method" would be a designation more 
nearly representing the facts of history than is "Newton's method." 

We proceed to the description of some details in Raphson's publica- 

* A copy of the edition of 1697 is in the New York Public Library. As this copy is in the reference depart- 
ment, I was not able to borrow it through a library loan. My knowledge of Raphson's book is derived mainly from 
the extracts given in the Latin edition of Wallis' Algebra, London, 1693, and the extracts given in an appendix, by 
Maseres, to W. Frend's Principles of Algebra, London, 1796, pp. 467, 468. 
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fir) 
tion. We have pointed out that the form now so familiar, i,\ \ , was not 

/ (r)' 

used by Newton, but was used by Raphson. To be sure, Raphson does not 
use this notation; he writes f(r) and /' (r) out in full as polynomials. Nor 
does Raphson find the first derivative by the rule of the calculus; his opera- 
tions are purely algebraical. He gives "canons" or "theorems" for the 
forms of the fractions, in which he writes down the polynomials /(r) and/' (r) 
for each equation up to that of the tenth degree inclusive. From the stand- 
point of algebraic notation it may be of interest to quote his expressions for 
one of these, say the quintic. Taking 1, b, c, d, ... as the coefficients of the 
given equation and g as the value of the approximation already reached, he 
writes for the quintic the following "canon:" 

"Pro potestate Quinta 

9 Q 9 g g 5 g g g g) 

bggg g Ubggg\ 
cggg 3cgg\x" 

d g g 2 dg 

fg f. 

The terms in the left hand column represent f{g) , those in the right hand 
column/' (g)x. With admirable consistency, Raphson gives the "canon," 
"pro potestate decima," by writing out each of the ten factors in gr 1 °. 

Newton and Raphson explained their methods only in connection with 
rational integral algebraic equations. The extension of the Newton-Raphson 
method to irrational and transcendental equations appears to have been 
made for the first time by Thomas Simpson, in his Essays... on Mathematics, 
London, 1740, p. 81. He does not mention Newton and Raphson, and calls 
his procedure a "new method." 

Nearly all eighteenth century writers and most of the early writers 
of the nineteenth century carefully discriminate between the method of 
Newton and that of Raphson. Then appear writers like Euler, Laplace, Lacroix 
and Legendre, who explain the Newton-Raphson process, but use no names. 
Finally, in a publication of Budan in 1807,* in those of Fourier of 1818 
and 1831, t in that of Dandelin in 1826, J the Newton-Raphson method is at- 
tributed to Newton. The immense popularity of Fourier's writings led to 
the universal adoption of the misnomer "Newton's method" for the 
Newton-Raphson process. 



* Nouvelle m&thode pour la resolution d. Squat, num., Paris, 1807, p. 78. 

t Bulletin d. sciences par la societe" philomatique de Paris, annee 1818, p. 61; Fourier, Analyse des equations 
dUerminies, 1831, pp. 169, 173. 177, etc. 

% N-. mSmoires d. I' acadimie r. d. sciences et belles-lettres de Bruxelles, T. 3, 1826, pp. 28, 29. 



